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This paper examines the consistency of the exact scaling laws for isotropic MHD tur- 
bulence in numerical simulations with large magnetic Prandtl numbers Pm and with 
Pm = 1. The exact laws are used to elucidate the structure of the magnetic and velocity 
fields. Despite the linear scaling of certain third-order correlation functions, the situation 
is not analogous to the case of Kolmogorov turbulence. The magnetic field is adequately 
described by a model of stripy (folded) field with direction reversals at the resistive scale. 
At currently available resolutions, the cascade of kinetic energy is short-circuited by the 
direct exchange of energy between the forcing-scale motions and the stripy magnetic 
fields. This nonlocal interaction is the defining feature of isotropic MHD turbulence. 



1. Introduction 

What is the structure of the saturated state of isotropic M HD turbulence? This is 
possibly the oldest question in the theory of MHD turbulence ( Batchelor 1950l ) — the 



answer to which is fundament ally impor tant to our understanding of cosmic magnetism 
(see, e.g, review bv [Schekochi hin fc Cow ley 2006,, and references therein). The problem 
can be posed in the following way. Consider the equations of incompressible MHD: 

du 

— +u-Vu = -Vp + B -VB + v/\u + (1.1) 
ot 

OB 

— + u VB = B Vu + rjAB, (1.2) 

ot 

where u is the velocity, B the magnetic field, p the total pressure determined by the in- 
compressibility condition V-tt = 0, / a body force, v viscosity and rj magnetic diffusivity 
(we use units in which p is sealed by p and B by ^/A^Tp, where p = const is the density 
of the medium). The evolution equations for the kinetic and magnetic energies are 



dt 2 

d{B^) 



-{BB ■.Vu)-iy{\Vu\^)+e, (1.3) 

2 {BB:Vu)-r^{\VB\^), (1.4) 

where (...) means volume averaging and e = (m • /) is the average injected power per 
unit volume (formally speaking, e depends on u and cannot be predicted in advance 
unless / is a white noise in time — in the latter case, e is fixed, which makes white-noise 
forcing an attractive modelling choice). The forcing acts at the system scale £o and we 
shall assume that the particular choice of / does not change the properties of turbulence 
at scales much smaller than £o. In hydrodynamics, one often considers decaying, rather 
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than forced, turbulence. The structure of the turbulence at small scales is expected to 
be the same for the decaying and forced cases. This appears not to be true for MHD, 
probably because of a tendency of the magnetic field t o decay towards large-sc ale force- 
free states (co mpare, e.g., the decaying si mulations o f Biska mp fc Miiller 2000l with the 
forced ones of ISchekochihin et aZ.ll2004l or lHaugen et a/,,, 2004) . We consider the case in 
which the velocity field is forced, while the magnetic field is not and can receive energy 
only via interaction with the velocity field. If there is dynamo action, an initially weak 
magnetic field is amplified by the turbulence until it becomes dynam ically significant an d 
saturates. Certain choices of the forcing, e.g., helical random forcing ( Brandenburg|[2001 ) . 
lead to the generation of magnetic field at scales larger than the forcing scale £o- In the 
presence of this mean field, the turbulence at the small scales is (globally) anisotropic 
— a distinct case that will not be discussed here (see review by ISchekochihin fc Cowlev 
Turbulence produced by a spatially homogeneous isotropic nonhelical random 



forcing does not generate a mean field, but, at least when the magnetic Prandtl number 
Pm = ly/v > 1, does lead to amplification of small-scale (r < ^o) magnetic energy (e.g.. 



Schekochihin et al\ l2004t iHaugen et al\ 120041 ) , a process known as small-scale dynamo. 



The saturated state of the small-scale dynamo is the fully developed isotropic MHD 
turbulence that is the subject of this paper. 

While some level of physical understanding of the small-scale dynamo and its satura- 
tion does exist ( ij3.ip . the detailed structure of the saturated state is still unknown. The 
existence of exact scaling laws analogous to Kolmogorov's 4/5 and Yaglom's 4/3 laws 
(equations (|2.2p and (|2.4p below) has occasionally been interpreted as suggesting that 
MHD turbulence has incrtial-range scaling analogous to the hydrodynamic case: specif- 
ically, that the fundamental physical fields in MHD are Elsasser variables ~ u± B, 
which have Kolmogorov scaling (from equation (|2.4p . the increments Sz^{r) ~ r^^^ for 
point separations r greater than the viscous and resistive scales), and also that the mag- 
netic and kinetic energies arc in scalc-by-scale equipartition, 6u(r) ^ 6B{r). Obviously, 
the exact laws by themselves do not rigorously imply any of this. The qualitative anal- 
ogy with Kolmogorov turbulence is based on the assumption that interactions occur 
between comparable scales (locality in scale space) and that the exact laws therefore 
describe a scale-by-scale energy cascade. However, in isotropic MHD turbulence, the in- 
teractions are not local. An obvious example of a nonlocal process is the small-scale dy- 
namo, whereby magnetic fields with direction reversals at the resistive scale (the folded 
fields) are generated by the velocity fluctuations that have much larger scales (i ^B.ip . 
Both numerical evidence and physical reasoning suggest that the saturated state of this 
process is controlled by the nonlocal interaction between the large- scale (forcing-scale) 



velocity gradients and the folded (small-scale) magnetic structures ijSchekochihin et 



[2004; A lexakis et al . 2005) • It is then an interesting question how such a saturated state 
can be consistent with the exact scalings mandated by the exact laws. Since these scal- 
ings represent the only rigorously established constraints on the statistics of the isotropic 
MHD turbulence, understanding how they are satisfied is a way to learn something about 
the structure of the turbulence. 

The purpose of this paper is to give an interpretation of the exact laws in various 
scale ranges, establish that the laws hold in numerically simulated MHD turbulence, and 
determine the extent to which the simulations access the desired asymptotic regimes. The 
plan of further developments is as follows. The derivation of the exact laws is reviewed 
in Sj2l A theoretical discussion of the generation of stripy (folded) magnetic fields and 
the structure of the turbulence at subviscous scales is given in fJS] Numerical tests are 
reported in 21 The concluding section, ^ discusses isotropic MHD turbulence in the 
inertial range and the unresolved issues. 
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2. Exact scaling laws for isotropic MHD turbulence 

The procedure for obtaining exact scaling laws is standard. Consider two points, x and 
x'. Denote it = u{x), u' = u{x'), Su = u'—u. Use equation (jl.ip taken at points x and x' 
to derive an evolution equation for the correlation tensor (uiu'j). This equation contains 
third-order correlation tensors {uiUju'/^) and {BiBju'j^). Because of the assumed spatial 
homogeneity, all two-point correlation tensors depend only on r = a;' — a;. We denote the 
projections of u and B on r by subscript l (for "longitudinal"), e.g., 6ul = Su-f, where 
r = r/r. Assuming isotropy, all correlation tensors can be expressed in terms of a set of 
scalar functions that depend only on r = \r \. The evolution equa tion of the second-order 



correlation function of the velocity field is (jChandrasekharjl 195 ih 



This is the von Karman-Howarth equation for MHD turbulence. At r = it reduces 
to equation (|1.3[) for the evolution of the kinetic energy. For r > 0, it contains infor- 
mation about the scale-by-scale energy budget (the turbulent cascade, energy exchanges 
between velocity and magnetic fields, etc.). If we now consider a statistically stationary 
state, in which all averages are constant in time, equation (j2.1ll transforms into a direct 
generalisation of Komogorov's 4/5 law for MHD turbulence: 

/ Q X / o V d(Suj) 4 

{Sul)-6{Bl6uL)-6i^^gf^^--er, (2.2) 

where we used the Taylor expansion {Jlu'j) = (l/3)e + 0{r'^ /l"^), in which only the first 
term needs to be retained for r £o- Note that the term involving mixed correlations 
of the velocity and magnetic fields contains the field B rather than its increment 5B — 
this hints at the nonlocality of the interaction between it and B, which, as explained in 

21 is the fundamental feature of MHD turbulence . 

Equ ation \2.2\ was derived by IChandrasekhar ( 195 if ) (see also Politano fc Pouquet 



199861 ). He also found the two other exact laws of MHD turbulence by constructing 
evolution equations analogous to equation (|2.ip for (B^i?^) and luiB'^') (representing the 
magnetic-energy and the cross- helicity budgets, respectively). These laws relate certain 
third-order mixed correlation functions of it and B to each other and to the viscous and 
resistive dissipation. The laws can be written in several equivalent forms, none of which 
is physically illuminating in any obvious way. It is, therefore, as good a choice as any to 
cast them in a form that directly generalises results known in fluid turbulence in the way 
equation (|2.2p does. This is achieved if we recall the alternative form of the equations of 
incompressible MHD in terms of the Elsasser variables = u ± B: 

^+z^-Vz^ = -Vp+^^Az^ + ^^Az^ + f. Vz^=0. (2.3) 
ot 



The form of the nonlincarity suggests a formal analogy with "passive" advection of the 
field z^ by the field z^ and vice versa, and it is, therefore, not surprising that developing 
evolution equations for the correlators (^z^z^'") leads to a pair of exact law s that are the 



MHD version of Yaglom's 4/3 law for passive scalar ( Politano fc Pouquet 1998ah 



{Sz^\Sz^\') - |: [(^ + V) {\Sz^\') + (^ - ^) {Sz+ ■ Sz-)] = -1 er, (2.4) 
where 6z^ = z^{x + r) — z^{x). 
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3. Theoretical considerations: linear stretching and stripy fields 

3.1. Small-scale dynamo and its saturated state 
The small-scale dynamo is due to random s tretching of the magnetic field by the ve- 



locity gradients asso ciated with turbulence jBatchelor 1950t Moffatt fc Saflman 1964 



Kazantsev Il968l : see ISchekochihin fc CowlevI l2006l for a review). In Kolmogorov tur- 



bulence, the velocity increments obey 5u{r) ~ (er)^/'^ for ^ ^ ^u, where £^ 
(z/3/e)'/^ ~ 4Re-3/^ is the viscous scale, Re = Wrms^o/^ is the Reynolds number and 
(u^Y^^ ~ (eeoy/^. Clearly, the viscous-scale motions give rise to the fastest 
stretching, with the rate 7 ^ 5u{ii,)/iy ~ (e/i^)^/^ ~ (Mrms/^o)R'e^^^. The random 
stretching produces magnetic fields folded in long flux sheets with fold length f y ~ In 
and direction reversals at the resistive scale ~ {vIlY^'^ ^ (.yFmT^^^ . When Pm 3> 1, 
the resistive scale is much smaller than the viscous scale i^. 

While the folded fields are formally small-scale, they are capable of exerting a back 
reaction on the flow that has long-scale spatial coherence — this is because the ten- 
sion force in equation (|1.2p is quadratic in B, so the direction reversals do not matter 
(note also that B ■ VB ^ independent of £^). iBatchelor ( 1950[ ) proposed that 



the small-scale dynamo would saturate when the magnetic field is strong enough to op- 
pose the stretching ac tion of the viscous-scale motions, ( B"^) ~ ieiy)^^'^ ^ (u^) Re~^^^ 
(see also Moffatt 1963 ) . The alternative view expressed by ISchliiter fc Biermann lll95d) 



was that the stretching action of the inertial-range motio ns would also be suppressed 
and (i?^) ^ (u^) in saturation. Schekochihin et al. ( 20041 ) argued that the latter view 



is more consistent with numerical evidence but that, contrary to Schliiter and Bier- 
mann's assumption of scalc-by-scale cquipartition. the magnetic field in the saturated 
state remained small-scale-dominated and folded, with folds elongating to the forcing 
scale, €|| ~ £0. They proposed that saturation is controlled by the balance of the stretch- 
ing of the field by the forcing-scale motions and the back reaction from the folded 
fields. The resistive scale is set by the characteristic time of the forcing-scale motions, 
£ri ~ {v^o/urms)^^^ ~ ^oRm"^^^. This can still be much smaller than the viscous scale 1,^. 
Let us now examine equations (|2.2p and (|2.4p in the subviscous scale range, r ^ £,j. 

3.2. The 4/5 law at the subviscous scales 

At subviscous scales, the velocity field is smooth, i.e., the velocity increments can be 
approximated by Taylor expansion: (5tt ~ r • Vtt(a;), where Vm does not depend on r. 
Let us substitute this into equation (|2.2p . The term (Suj^) ^ is subdominant. Since 

(Sul) ~ ((rr : Vtt)') = {V,UjVkUi) {hfjhn) = ^ (I Vwp) r^, (3.1) 

the viscous term is (4/5)i^ (| Vitp) r. The magnetic term is 

6 (bISul) 6 {BkBiV^Uj) {hfjhri) r = \ {BB : Vu) r. (3.2) 

In equations ((3Tl) - ((3^ . we used {hfjfkh) = {SijSki + SikSji + SuSjk). We see 

that equation p.2p reduces to the power balance 

e= (SB : Vm) +zy(|VMp), (3.3) 

which is the steady-state form of equation (jl.3p . Thus, the 4/5 law carries no new infor- 
mation — it simply registers the fact that the injected kinetic energy is partly dissipated 
viscously, partly transferred into the magnetic field. The latter part is dissipated resis- 
tively: {BB : Vit) = r/(|VSp) from the steady-state form of equation (|1.4p . 
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Figure 1. (a) The stripy field model. (&) A cross section of for the run used in figure [JJ a) 
(Pm = 1250). Arrows indicate the in-plane direction of the field. 



3.3. The 4/3 law and the stripy-field model 
The 4/3 law, equation (|2.4[) . is more interesting. In terms of Su and SB, it reads 

+ {5ul\SB\'') - 2 (JSi^M • SB) T {SBlISB^) 
= er + 2^ [^{\Su\'-) ± + 77) (<5« • SB) + v{\SB\^)] . (3.4) 

Since Su ^ r, the first term on the left-hand side is r^, the second and third terms are 
~ r^, so all three are subdominant and can be dropped. The remaining terms depend on 
the structure of the magnetic field. Let us consider a drastically simplified model of this 
structure based on the understanding that the field is organised in folds with direction 
reversals at the resistive scale. The model is illustrated in figure [TJ a). The orientation of 
B on this sketch must be understood as one among many possible orientations with equal 
probability, since the flow is statistically isotropic. The field is taken to be straight (or 
coherent along itself on scales comparable to the forcing scale £0) and to oscillate rapidly 
with spatial period 2^^ <C io across itself. Given an arbitrary point x, if we randomly 
pick a point x' on a sphere centered at x with a fixed radius r such that £jj r <^ i^, the 
field increment will be SB = —2B{x) with probability 1/2 and SB = with probability 
1/2. We may now use this simple model and {rifj) = {l/3)Sij to calculate 

{Sul\SB\'') =2{B^V,Uj) ihf,) r = ^ {B^V • «) r = 0, (3.5) 

2 {SBlSu -SB) ^4 {B.BjVkUj) {hfk) r = ^ {BB : Vu) r, (3.6) 
{SBl\SB\^) = -4{B'B,) in) = 0, (3.7) 

(|(5n|2) = {W,u,WkU,) {hh) r^ = l {\Vu\^) r\ (3.8) 
{Su ■SB) = - {B,V,u,) (fj) r = 0, (3.9) 
{\SB\^) = 2{B^). (3.10) 
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Substituting these expressions into equation ()3.4p . we see that it reduces to the power 
balance p.3|) just as the 4/5 law did. It is essential, however, that to achieve this seemingly 
trivial outcome, we had to make assumptions about the spatial structure of the magnetic 
field. We have thus shown that a stripy alternating field is consistent with the 4/3 law. 

It is perhaps appropriate to stress that the stripy field sketched in figure[lja) is a highly 
idealised model of the field structure that does not take into account, for example, such 
features as turning points, bending of the folds on the scale of the flow, the extent to which 
the stripes are volume filling, possibility of some degree of misalignment of the alternating 
fields, etc. However, it has allowed us to give a particularly simple demonstration of the 
way the 4/3 law accommodates the folded magnetic fields evident in figure [ij 6). The 
key to understanding why such a simple model appears to be sufficient at least on the 
qualitative level may be that the fields best described by this model are also the strongest 
ones, so the model captures well the statistical quantities weighed towards the parts of 
the system where the field strength is largest. 



4. Numerical Tests 

We now consider how equations (j2.2p and (j2.4p are satisfied in numerical experiments. 
We use the data f r om d irect numerical simulations of isotropic MHD turbulence by 



Schekochihin et al\ (|2004l ). The list of runs and other details are provided in that paper. 
These are three-dimensional spectral simulations in a periodic cube of size 1. They employ 
a white-noise random forcing at the box scale. The average power input is e = 1. 

4.1. The averaging procedure 

We use three levels of averaging to obtain two-point correlation and structure functions. 

• Spherical average. For a given point x and a given point separation r, we average 
over all possible orientations of r. We use 200 discrete orientations. In order to obtain 
accurate results, it is necessary to resort to linear interpolation between grid points to 
compute the fields at x' = x + r. Using the spherical average ensures that spurious 
anisotropics introduced by finite size Cartesian grid effects are filtered out. It is not 
restrictive, because the turbulence is statistically isotropic. 

• Volume average. Since the turbulence is spatially homogeneous, we average over all 
points X. In practice, the points x are picked on a submesh that is 10 times sparser than 
the mesh of the simulations. 

• Time average. Finally, we average over about 20 box-crossing times. This is equiv- 
alent to averaging over many realisations of the turbulence. A relatively large number 
of independent snapsh ots is necessa ry in order to obtain converged correlation func- 
tions of odd power (see Rincon 2006h . The averaged quantities we report arc sufficiently 



converged (within a few percent for third-order quantities, which have the slowest con- 
vergence rate) for a qualitative understanding of the scale-by-scale budgets. We do not 
attempt any high-precision determination of scaling exponents from this data. 

Note that in computing the two-point statistical quantities, we may only use points 
x' within a radius r up to 1/4 of the box size from any given point x because at larger 
distances spurious correlations arise due to the periodicity of the box. 

4.2. Numerical results in the viscosity- dominated limit 

It is not currently possible to have a resolved numerical simulation with Pm ^ 1 and 
Re 2> 1. However, one can study the subviscous fields in numerical simulation s with 
Pm ;:g> 1 and Re ~ 1 (the viscosity-dominated limit; see ISchekochihin et al. 2004 ). This 



means that the inertial range physics is sacrificed and we study the saturated state of a 
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Figure 2. (a) The 4/3 law for a run with Pm — 1250, Re ~ 1 (run S5 oflSchckoch ihin et aH20M : 
256^, f = 5x 10"^, Tj = 4x10"^^). We plot minus the left- (+++) and the right-hand ( — ) sides of 
equation (|2.4p , and also minus the dissipative term (•■■). Several individual third-order structure 
functions in the left-hand side of equation (|3.4|l are shown: — (^Sz'[\dz\^'j ), 2 (SBlSB ■ Su) 

{— ■ — ), — ^(5ii£ ( ), — (^Sux^lSBl^'j (--■■■ — ). The three functions that are not shown 

are small. (6) The second-order structure function (\SB\'^) /2 (B^) for a seq uence of runs with 
increasing Pm and constant Re ~ 1 (runs S1-S6 of ISchekochihin et a/.|[2003 ). 



small-scale dynamo with Rm ^ 1 and a singl e-scale (^ £n) ran doml y forced flow. We note 



that in this setting, the difference between thc'Ba tcheloi (1950) and Schliiter fc Biermann 



( 19501) scenarios of saturation mentioned in > j3.1l cannot be detected 



Since the viscosity is large, the flow is smooth at all scales. As we explained in 
the 4/5 law in this regime does not contain any information beyond the power bal- 
ance p.3p . Numerical results are in line with this expectation (not shown). The con- 
sistency of the 4/3 law is illustrated in figure [21 a), which shows that the left- and 
right-hand sides of equation (j2.4[) independently computed from the numerical data 
match each other. The deviations at large scales are due to the higher-order terms in 
the Taylor expansion (/ • u') ~ e + 0{r'^ /£q). In equation (j2.4p . the dominant term 
in the left-hand side is the viscous one. However, for a smooth velocity held, it is lin- 
ear in r (equation p.Sp ). and, indeed, wc find that the third-order structure function 
— (Jz^jiJ^+p) ~ (4/3) (e — j/(|Vitp)) r is also linear in r for r > The main con- 
tribution to it is due to 2 {SBlSu ■ SB), while all other terms in the left-hand side of 
equation p.4p are subdominant. This is consistent with what is predicted by the stripy- 
field model, equations p.5p - (|3.7p . although equation (|3.6p is only satisfied approximately 
(within ~ 30%). Note that for r £jj, the magnetic field is smooth, SB ^ r, so the tran- 
sition from the subresistive scaling (5z£|(52;+p) ~ r^ to <^(5z£|(52;+p) ~ r shows what the 
effective magnitude of the resistive scale is: 4; ~ 0.03 for the run used in figure [2j a). 

4.3. Scaling of the magnetic-field increments and the magnetic- energy spectrum 

The stripy-field model implies that 6B{r) ^ i?rms = const for <C r ^ £^ (equa- 
tion ((3T0)) ). Note that {\SB\'^) ~> 2{B'^) in the limit r oo for any field structure, 
as long as the field decorrelates at long distances, (S' • B) — > 0. The nontrivial prop- 
erty of the stripy-field model is that this value is reached already at r > £^ because 
(B' • B) = due to direction reversals. Figure [DJ 6) shows that the scale at which the 
value (|i5Bp) ~ 2 (i?^) is reached does indeed decrease as Rm is increased. 

What the magnetic-energy spectrum is in the saturated state is not currently known. 
The stripy-field result SB ~ r° appears to imply that the spectrum is k~^. However, the 
stripy-field model does not take into account that the alternating fields do not cancel 
each other perfectly, so there should be a gradual loss of correlation as r increases. 
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Figu re 3. (a) The 4/5 law for a run with Pm = 1, Re ~ 400 (run A of ISchekochihin et al\ 
120041 : 256^, = rj ^ 5 x 10""). We plot minus the left- (+++) and the right-hand ( - ) sides of 

equation (|2.2|l . minus the viscous term (• ■ ■), — ( ) and 6 (^B\5ul) (— ■ — )• Inset: plot 

of £t = —6 (^B\5ul) r I (^Su^). (b) The 4/3 law for the same run. See caption of figure [2j a) for 
explanation of line styles. The functions {SulSu ■ SB), (^6Bl\Su\'^) and (SBl\SB\'^) are again 
small (not shown). 



If this effect is important, the correlation function may have a peak around 

r ^ £ri followed by a gradual fall off to the value 2 (i?^) at r ^ If this fall off is a 
power law ~ {^■q/i")^'^" , the spectrum should then be fc". For example, in the kinematic 
regime of the small-scale dynamo, theory predicts a = 3/2 |Kazantsev|[l968i) . In the 
saturated state, the numerically computed correlation functions (figure [2^6)) show no 
sign of having a peak, which s uggests a = —1. It cannot, however, be ruled out that, 
as suggested in ISchekochihin et a l. (200^. the spectrum is flatter, possibly even a > 0, 
but the convergence to the asymptotic scalings is slow and the values of Pm used in our 
simulations are insufficient to determine a accurately. 

4.4. Numerical results for Pm = 1 

While the limit Rm Re ^ 1 is unresolvable, some idea of the structure of isotropic 
MHD turbulence in a regime that is not viscosity-dominated can be gained by studying 
the case of Pm = 1 , which has been the favourite choice in numerical experiments (e.g., 
Haugen et al Such simulations may, in fact, be more relevant than it would appear 



to understanding the Pm ^ 1 regime if the saturated state of the small-scale dynamo is 
controlled by the interac tion of the f olded magnetic fields with the forcing-scale velocities 
( ^3.ip . rather than — as Batchelor ( 1950f ) thought — with the viscous-scale ones. 



Whether this is true is best addressed numerically by comparing the four terms in 
the 4/5 law, equation (j2.2p . This is done in figure [3lJ a) for a run with Re = Rm ~ 400. 
The 4/5 law is well satisfied: the independently computed left- and right-hand sides 
of equation (|2.2p match except for the large-scale discrepancy due to the higher-order 
contributions from the forcing term. The viscous term is subdominant for r > ~ 0.05 
leaving some room for a possible inertial range. Except at r ~ the dominant balance 
is between Q(^B\5ui,) and (4/5)er, while {5u\) is subdominant. Physically this means 
that the kinetic- energy cascade is short-circuited by the energy transfer into the magnetic 
field. Interestingly, for the subdominant "cascade" term, we find empirically 

^{5ul)^^Bl5u,)L^\^r\ (4.1) 

where ~ 0.25 is some typical length, comparable to the forcing scale (.q. The resolution 
is limited and the scale range where {Suj^^ ex holds is not wide. Numerically, the 
relation oc {B\5uL)r turns out to be satisfied across a somewhat wider interval 
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(see inset in figure [3^ a)). Equation (|4.ip suggests that the velocity increments have the 
scaHng Su ^ r^/"^. Indeed, the second-order structure function of it computed from our 
data is consistent with {5u\') ^ r'^l'^ (not shown) and tlie kinetic- energy spectrum does 
appear to have the corresponding scahng A:"^/"^, as noticed bv lSchelcochihin et al\ (|2004h . 

The 4/3 law is well satisfied (figure [3^6)), as it was in the viscosity-dominated case 
(i j4.2p . but now there is an interval (r > 0.02) where the dissipative terms are small and 
equation (|2.4p becomes — {8z'^\8z^^') ~ (4/3)er. The numerically computed third-order 
structure function in the right-hand side of this equation does, indeed, scale linearly 
with r. The main contribution to this function still comes from 2 {SBlSu ■ SB) (see 
equation p.4p ). except at the largest scales, where (^Sul\Su\'^') becomes important. That 
the 4/3 law approximately reduces to a balance between the magnetic-interaction term 
2 {SBlSu ■ SB) and (4/3)er with the kinetic-energy "cascade" term interfering 
only at r ~ £o is in line with what we learned above from the 4/5 law. The qualitative 
similarity of the internal structure of {5z2\Sz'^\'^') (the relative importance of the six 
terms in the left-hand side of equation p.4p ) to the viscosity-dominated case appears to 
suggest that the magnetic field in the Pm = 1 case still has a stripy (folded) stru cture. 
This is, indeed, supported by various quantitative tests (jSchekochihin et aZ.1 120041 ). 



5. Discussion: turbulence in the presence of stripy field 

Numerical results appear to make a fairly compelling case that the structure of the 
magnetic field in the isotropic MHD turbulence is dominated by folds (stripes). We have 
demonstrated above that such a structure is consistent with the exact scaling laws. It has 
proven to be much more difficult to understand the structure of the velocity field. It is 
clear that the cascade of kinetic energy is short-circuited at the forcing scale, with injected 
energy diverted into maintaining the folds against continuous Ohmic dissipation. It is not 
clear from the available numerical results whether most or only a finite fraction of the 
injected power e is thus diverted. If the latter is the case, i.e., if the magnetic-interaction 
term & (^B\5ui,) only cancels a part of the flux term (4/5)er in equation (|2.2p . there 
should still be a kinetic-energy cascade from the forcing to viscous scale and it is hard to 
see from equation (|2.2p how the velocity increments associated with this cascade can fail 
to have the Kolmogorov scaling, {Su\) ~ — (4/5)xe?', where x < 1 is some finite positive 
number. These motions a re likely to be Alfvenic perturbations of the folded structure 
( Schekochihin et al ] |2002|) . If, on the other hand, the short-circuiting of the kinetic-energy 
cascade is nearly complete, i.e., Ql^B^Sui,) ~ (4/5)er, then a subdominant scaling is 
possible: {Su^^ ~ r^, /3 > 1. Equation ()4.ip . which we deduced from our numerical 
experiment, appears to support this possibility with /3 = 2. However, the resolution of 
our study was not sufficient to have a convincing parameter scan in Re and prove that 
the steeper-than-Kolmogo rov scaling of the vel ocity increments is not simply due to too 
much viscosity. Note that iHaugen et al\ (j2004l ) report kinetic- energy spectra that may 
be consistent with k~^/^ at values of Re roughly twice as large as ours. 

We remark in passing that the short-circuiting of the kinetic- energy cascade appears 



also to be a fea ture of the turbulence in polymer solutions (jPe Angelis et all 12005 



iBerti et~al. 20061 ). the mathematical description of which is similar to MHD, with elastic 
polymer chains stretched by turbulence in the same way magnetic fields are. 

While the above study does not constitute the complete solution of the problem of 
isotropic MHD turbulence, it adds to the weight of evidence that this turbulence is con- 
trolled by the direct nonlocal interaction between the forcing-scale motions and small-scale 
stripy, or folded, magnetic structures. This nonlocality means that scaling arguments 
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cannot be made along the usual, Kolmogorov-like, lines. This makes the problem con- 
ceptually difhcult. The numerical simulations, while helpful, are not as yet sufficiently 
large to access the asymptotic regime of interest: ideally, Rm ^ Re ^ 1. The exact 
scaling laws provide the only rigorous quantitative constraints available in the theory of 
turbulence and arc, therefore, useful. In this paper, we have attempted to discern the 
message that these laws carry about the structure of isotropic MHD turbulence. Given 
the limited resolution of the numerical experiments that we used to guide us, it may be 
fruitful to return to this type of analysis when larger computations become feasible. 
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